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Abstract 

The doping-dependent evolution of the d-wave superconducting state is studied from the perspective of the 
angle-resolved photoemission spectra of a high-7c cuprate, Bi2Sr2CaCu208+5 (Bi2212).The anisotropic evolution of 
the energy gap for Bogoliubov quasiparticles is parametrized by critical temperature and superfluid density. The 
renormalization of nodal quasiparticles is evaluated in terms of mass enhancement spectra. These quantities shed 
light on the strong coupling nature of electron pairing and the impact of forward elastic or inelastic scatterings. We 
suggest that the quasiparticle excitations in the superconducting cuprates are profoundly affected by 
doping-dependent screening. 
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Background 

Electronic excitations dressed by the interaction with the 
medium are called quasiparticles. They serve as a direct 
probe of the anisotropic order parameter of a super- 
conducting phase and also as a clue to the electron- 
pairing glue responsible for the superconductivity. In 
fact, the major unresolved issues on the mechanism of 
high-Tc superconductivity depend on the low-energy 
quasiparticle excitations. The superconducting order 
parameter, which is typified by the particle-hole mix- 
ing and gives rise to Bogoliubov quasiparticles (BQPs), 
manifests itself as an energy gap in quasiparticle excita- 
tion spectra. In cuprate superconductors, however, the 
energy gap increases against the decrease in critical tem- 
perature Tc with underdoping and is open even at some 
temperatures above Tc [1-3]. In the direction where the 
(i-wave order parameter disappears, renormalization fea- 
tures have been extracted quantitatively from the gapless 
continuous dispersion of nodal quasiparticles (NQPs), 
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suggesting strong coupling with some collective modes 
[4]. Nevertheless, the origins of these features remain 
controversial [4,5]. 

In this paper, we address the doping dependence of 
BQP and NQP of a high-Tc cuprate superconductor, 
Bi2Sr2CaCu208+5 (Bi2212), on the basis of our recent 
angle-resolved photoemission (ARPES) data [6-8]. The 
use of low-energy synchrotron radiation brought about 
improvement in energy and momentum resolution and 
allowed us to optimize the excitation photon energy. After 
a brief description of BQP and NQP spectral functions, 
we survey the superconducting gap anisotropy on BQPs 
and the renormalization features in NQPs. In light of 
them, we discuss possible effects of doping-dependent 
electronic screening on the BQP, NQP, and high-Tc 
superconductivity. 

Methods 

High-quality single crystals of Bi2212 were prepared by 
a traveling- solvent floating-zone method, and hole con- 
centration was regulated by a post-annealing procedure. 
In this paper, the samples are labeled by the Tc value in 
kelvin, together with the doping-level prefix, i.e. under- 
doped (UD), optimally doped (OP), or overdoped (OD). 
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ARPES experiments were performed at HiSOR BL9A in 
Hiroshima Synchrotron Radiation Center. The ARPES 
data presented here were taken with excitation-photon 
energies of hv = 8.5 and 8.1 eV for the BQP and NQP 
studies, respectively, and at a low temperature of T = 
9 — 10 K in the superconducting state. Further details of 
the experiments have been described elsewhere [7-9]. 

The relation between a bare electron and a renormal- 
ized quasiparticle is described in terms of self-energy 
^k(0> which can be regarded as a factor of feedback on 
the wave function from past to present through the sur- 
rounding medium. Incorporating a feedback term into the 
Schrodinger equation, we obtain 



dt 



J -oo 



(1) 



where V^k(0 ^rid co^ denote a wave function and a bare- 
electron energy, respectively. It is obvious from Equation 1 
that the self-energy is a linear response function. There- 
fore, its frequency representation, ^i^ico), obeys the 
Kramers -Kronig relation. As the solution of Equation 1, 
we obtain the form of dressed Greens function. 



(2) 



The spectral function given by Ai^(cjo) = — Im G\^{oj)/7Z is 
directly observed by ARPES experiments. The extensive 
treatments of the ARPES data in terms of Green's function 
are given elsewhere [10]. 

Results 

Superconducting gap anisotropy 

In the superconducting state, the condensate of electron 
pairs allows the particle-like and hole-like excitations 
to turn into each other. Hence, the wave function 
of a hole-like excitation also comes back by way of 
the particle-like excitation. Such a feedback has a 
sign-reversed eigenenergy, —oj^y and is expressed by 

EkCO = -/0(O A^e^'^^e-^^ - /Tk^a), where ©a), 
Ak, }\ and denote the step function, the particle-hole 
off-diagonal element, and the scattering rates of the inter- 
mediate and bare-particle states, respectively. The Fourier 
transform of gives the frequency representation of 
the self-energy of the BQPs, 



(3) 



Figure 1 shows the ARPES spectra of BQPs for under- 
doped and overdoped Bi2212 samples with = 66 and 
80 K (UD66 and OD80, respectively) [8]. As shown in 
Figure lb,c, an energy distribution curve was extracted 



from the minimum gap locus for each off-node angle 0 
and symmetrized with respect to the Fermi energy co = 0. 
These spectra were well fitted with a phenomenological 
function. 



(4) 



except for a featureless background. Equation 4 is deduced 
from Equation 3 and co^ = 0, neglecting yk after Norman 
et al. [11]. Figure lb,c exemplifies that the superconduct- 
ing gap energy A at each 0 is definitely determined by 
sharp spectral peaks. In Figure ld,e, the obtained gap 
energies (small yellow circles) are plotted over the image 
of spectral intensity as a function of sin 2^, so that the 
deviation from a <i-wave gap is readily seen with refer- 
ence to a straight line. While the superconducting gap 
of the overdoped sample almost follows the /i-wave line, 
that of the underdoped sample is deeply curved against 
sin 20, Furthermore, Figure Id indicates that the deviation 
from the (i-wave gap penetrates into the close vicinity of 
the node and that it is difficult to define the pure <i-wave 
region near the node. Therefore, the next-order harmonic 
term, sin 6^, has been introduced, so that the smooth 
experimental gap profile is properly parametrized [12-14]. 
The next-order high-harmonic function is also expressed 
as A(0) = An sin 2^ + (A* - AN)(3sin26> - sin66')/4, 
where the antinodal and nodal gap energies are defined as 
A* = A(^) 15/^450 and An = ^(dA/dO)\^^Q, respectively, 
so that An/ A* = 1 is satisfied for a pure <i-wave gap. 

The doping dependences of the superconducting gap 
parameters are summarized in Figure 2. One can see 
from Figure 2a that as hole concentration decreases with 
underdoping, the nodal gap energy 2 An closely follows 
the downward curve of 8,5kBTc in contrast to the mono- 
tonic increase in the antinodal gap energy 2A*. It seems 
reasonable that Tq primarily depends on An rather than 
A* for the underdoped Bi2212, because it follows from 
2 A* ^ 4ikBTc that the thermal quasiparticle excitations 
concentrate in the vicinity of the node and hardly occur 
around the antinode. The relevance of the nodal exci- 
tations has also been suggested by various experiments 
[15-19]. Then, the problem with Tc is that the nodal gap 
An is suppressed relative to the antinodal gap A*. This 
behavior can be associated with low superfluid density Ps 
[20]. Figure 2b,c shows that the doping dependence of the 
nodal-to-antinodal gap ratio An/ A* is quite similar to 
that of the square-root superfluid density [8,21,22]. 
The normalized gap plot in Figure 2d indicates that what 
occurs with underdoping is analogous to the nodal gap 
suppression observed with increasing temperature [17] 
in terms of the decrease in ps- It is notable that the 
square-root dependence on ps is a typical behavior of the 
order parameter as expected from the Ginzburg-Landau 
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Figure 1 Superconducting gap manifested in BQP spectra. The data are for underdoped and overdoped Bi221 2 samples with 7c = 66 and 80 K 

(UD66 and OD80, respectively) [8]. (a) Momentum-space diagram for an off-node angle, 0, and a bonding-band (BB) Fermi surface along which the 
ARPES spectra were taken, (b, c) Symmetrized energy distribution curves (colored circles) and their fits (black curves), (d, e) ARPES spectral images 
as a function of energy co and sin 26. Superimposed are the gap energies (yellow circles) and high-harmonic fit (yellow curve) as functions of sin 26. 



theory [23]. These findings can be written down in a 
simple relational formula, 

8.5/:Brc = 2AN = A*/5' (5) 
V Ps 

where :^ 31 |am~^, for a wide hole-concentration range 
ofBi2212. 

As presented in Figure 2e, the correlation between 
the nodal and antinodal gaps provides a perspective 
of crossover for our empirical formula (Equation 5). 
It is deduced from the conventional Bardeen-Cooper- 
Schrieffer (ECS) theory that l/^/k^T^ = 4.3 in the weak 
coupling limit for the <i-wave superconducting gap [23]. 
However, the critical temperature Tc is often lower than 
that expected from the weak coupling constant and a 
given A as an effect of strong coupling. Thus, the gap-to- 
Tc ratio is widely regarded as an indicator for the coupling 
strength of electron pairing and adopted for the coordi- 
nate axes in Figure 2e. As hole concentration decreases 



from overdoped to underdoped Bi2212, the experimen- 
tal data point moves apart from the weak coupling point 
toward the strong coupling side, and a crossover occurs 
at 8.5, which is about twice the weak coupling constant. 
It appears that the evolution of An is confined by two 
lines as An < 0.87 A* and 2 An < 8.5/cBrc- As illustrated 
in the insets of Figure 2e, the strong coupling allows the 
electrons to remain paired with incoherent excitations. As 
a result, the superconducting order parameter is reduced 
with respect to the pairing energy. Indeed, it has been 
shown that the reduction factor due to the incoherent pair 
excitations has a simple theoretical expression ^ps/pf^^ 
and that the nodal and antinodal spectra are peaked at 
the order parameter and at the pairing energy, respec- 
tively, taking into account a realistic lifetime effect [24,25]. 
Therefore, the latter part of Equation 5 is consistent with 
the strong coupling scenario, and furthermore, the two 
distinct lines in Figure 2e are naturally interpreted as the 
energies of the condensation and formation of the electron 
pairs. 
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Figure 2 Doping dependences of superconducting gap parameters, (a) Nodal gap energy 2An (blue circles) and antinodal gap energy 2A* 
(red squares) [8]. The solid curve denotes an energy of 8.5/cb7"c. (b) Square of nodal-to-antinodal gap ratio (An/A*)^ determined from ARPES [8]. 
(c) Superfluid density ps determined from magnetic penetration depth (triangles) [21] and from heat capacity (crosses) [22]. (d) Superconducting 
gap profiles normalized to the antinodal gap for underdoped and optimally doped samples with Tc = 42, 66, and 91 K (UD42, UD66, and 0P91 , 
respectively), (e) Correlation between 2An//cb7"c and 2A*//cb7"c ratios. The insets illustrate the occurrence of incoherent electron pairs in strong 
coupling superconductivity. 



Renormalization features in dispersion 

In the nodal direction where the order parameter disap- 
pears, one can investigate the fine renormalization fea- 
tures in dispersion. They reflect the intermediate-state 
energy in coupling between an electron and other exci- 
tations, and thus provide important clues to the pairing 
interaction. As for the electron-boson coupling, the inter- 
mediate state consists of a dressed electronic excitation 
and an additional bosonic excitation (Figure 3a). Averag- 
ing the momentum dependence for simplicity, the energy 
distribution of the intermediate state is expressed by 
A((jo - Q) S((jo -Q) -\- A(a) + Q) S(-(jo - Q) for a given 
boson energy Q and for zero temperature, owing to the 
Pauli exclusion principle. Therefore, taking into account 
the effective energy distribution of the coupled boson, 
a^F(Q), the self-energy is written down as follows: 

X \A(a)'-^)@(a)'-^) -\- A(a)'-\-^)&(-a)'-^)\ (6) 



J CO - co' -\- iO+ Jo 

X ^A(a)'-^)@(a)'-^) -\- A(a)'-\-^)G(-a)'-^)^, 

(7) 

where 0+ denotes a positive infinitesimal. 

In ARPES spectra, the real and imaginary parts of self- 
energy manifest themselves as the shift and width of 
spectral peak, respectively. Specifically, provided that the 
momentum dependence of ^i^(co) along the cut is negli- 
gible, and introducing bare electron velocity vq by co^ = 
vq /c, it follows from Equation 2 that the momentum distri- 
bution curve for a given quasiparticle energy co is peaked 
at k(co) = [co — ReE(a;)]/vo and has a natural half width 
ofAk(co) = -ImE(a;)/vo. 

We argue that the mass enhancement function defined 
as the energy derivative of the self-energy, X(co) = —(d/dco) 
is useful for the analysis of NQP [7,26]. The real 
and imaginary parts of k{co) are directly obtained from 
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Figure 3 Simulation for a single coupling mode at = 40 meV. Dotted and solid curves denote those with and without a d-wave gap of 

A = 30 meV, respectively, (a) Diagram of electron-boson interaction, (b) Eliashberg coupling function a^F(-(jo), dispersion kico) = [a;+ReE(a;)]/\/o, 

and momentum width Akico) = -lmE(a;)/\/o. (c) Real and imaginary parts of 1 + Xico). 



the ARPES data as the inverse of group velocity, Vg(co), 
and as the differential scattering rate, respectively. 

-^/cM = - [l + ReHoj)] = (8) 
d ^ I 

--Ak(co) = — Imk (co) (9) 
aco vq 

We note that —lmk(co) represents the energy distribution 
of the impact of coupling with other excitations and can 
be taken as a kind of coupling spectrum. However, it 
should be emphasized that —lmX(co) is expressed as a 
function of quasiparticle energy co, whereas the widely 
used Eliashberg coupling function a^F(Q) is expressed as 
a function of boson energy Q. For example, a simulation 
of k(co) using Equations 7 to 9 is presented in Figure 3b,c, 



where a single coupling mode is given at ^2 = 40 meV. 
One can see that the peak of a^F(—co) is reproduced 
by —lmk(co), provided that A(co) is gapless and approx- 
imated by a constant. As an energy gap of A opens in 
A(co), the peak in — ImA.(a;) is shifted from Q into ^2 + A. 
Nevertheless, irrespective oiA(co), the causality of ^(co) 
is inherited by k{co), so that Rek(co) and lmk{co) are 
mutually convertible through the Kramers-Kronig trans- 
form (KKT). The directness and causality of X(co) enable 
us to decompose the quasiparticle effective mass without 
tackling the integral inversion problem in Equation 7. 

Figure 4 shows the ARPES spectra along the nodal 
cut perpendicular to the Fermi surface for the super- 
conducting Bi2212 [7]. Although the splitting due to the 
Cu02 bilayer is minimum at the nodes, it has clearly 
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Figure 4 Dispersion Icinlcs manifested in NQP spectra. The ARPES spectra were taken in the superconducting state for Bi221 2 [7]. 

(a) Underdoped sample with 7c = 66 K(UD66). (b) Optimally doped sample with 7c = 91 K (0P91). (c) Overdoped sample with 7c = 80 K(OD80). 
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been observed by using some specific low-energy photons 
[6-8]. A prominent kinl< in the NQP dispersion is observed 
at 65 meV for all the doping level, as has been reported 
since early years [4]. In addition to this, another small kink 
at 15 meV is discernible in the raw spectral image of the 
underdoped sample (UD66) [7,27]. 

The fine renormalization features in the NQP disper- 
sion were determined by fitting the momentum distri- 
bution curves with double Lorentzian. Figure 5a,d shows 
the real and imaginary parts of k(co)/vo experimentally 
obtained as the energy derivatives of the peak position and 
width, respectively. The KKT of ReA.(a;)/vo in Figure 5a 
is shown in Figure 5b as ImA-(a;)/vo, which is comparable 
with the data in Figure 5d. A step-like mass enhancement 



in Figure 5a and a peak-like coupling weight in Figure 5b,d 
are consistently observed at 65 meV. This is a typical 
behavior of the mode coupling, as shown by the simula- 
tion in Figure 3. It is also found that an additional feature 
around 15 meV is dramatically enhanced with underdop- 
ing. In order to deduce the partial coupling constant, we 
express the mass enhancement factor A. as the form of 
KKT, 

X = ReA.(0) = - / —dco, (10) 

n Jo 0) 

Dividing the energy range of the integral in Equation 10, 
one can quantify the contribution from a particular energy 
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Figure 5 Doping dependences of NQP properties. The real and imaginary parts of mass enliancement spectra were directly deduced from the 
APRES data shown in Figure 4 [7]. (a) Inverse group velocity, 1 /Vgico) = [1 + Rek^co)] /vq, determined from (d/dco) k(co). (b) Differential scattering 
rate -\rf]X((jo)/vo, deduced from the Kramers-Kronig transform (KKT) of (a), (c) Partial coupling constants, X^^ (red circles) and A.'^ (blue triangles), 
deduced from (b). Also shown are the inverse group velocities at a; = 0 (black circles) and at a; = -40 meV (black triangles), (d) Differential 
scattering rate -Im X{co)/vq, directly determined from -(d/dco) Akico). (e) Partial coupling constants, X^^ (red circles), ^'^ (blue triangles), and X^^ 
(green diamonds), deduced from (d). 
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part. We refer to the KKT integrals of lmX(co)/vo for 
the low-energy (LE; 4 < < 40 meV), intermediate- 
energy (IE; 40 < \(o\ < 130 meV), and high-energy 
(HE; 130 < \co\ < 250 meV) parts as X^^/vq (red 
circles), X}^/vo (blue triangles), and k™/vo (green dia- 
monds), respectively. Those obtained from the data in 
Figure 5b,d are plotted in Figure 5c,e, respectively. Also 
shown in Figure 5c are the inverse group velocities at 
0) = 0 meV (black circles) and Btco = —40 meV (black tri- 
angles). Figure 5c and Figure 5e consistently indicate that 
as hole concentration decreases, the contribution of the 
low-energy part rapidly increases and becomes dominant 
over the other parts. 

Possible origins of the low-energy kink are considered 
from the energy of 15 meV and the evolution with under- 
doping. The quasiparticles that can be involved in the 
intermediate states are limited within the energy range 
of 1 60 1 < 15 meV, and the irrelevance of the antin- 
odal states is deduced from the simulation in Figure 3c. 
Therefore, the low-energy kink is due to the near-nodal 
scatterings with small momentum transfer. The candi- 
dates for bosonic forward scatterers are the low-frequency 
phonons, such as the acoustic phonons and the c-axis 
optical phonons involving heavy cations [7,28-31]. On the 
other hand, it has also been argued that the elastic forward 
scattering by off-plane impurities may give rise to the 
low-energy kink for the (i-wave superconductors [7,32]. 
In usual metal, both the potentials of the low-frequency 
phonons and the static impurities are strongly screened by 
the rapid response of electronic excitations. Therefore, the 
enhancement of the low-energy kink suggests the break- 
down of electronic screening at low hole concentrations 
[7,28]. 

The dispersion kink at 65 meV has been ascribed to an 
intermediate state consisting of an antinodal quasiparticle 
and the Big buckling phonon of ^2 ~ 35 meV [33]. How- 
ever, the mass enhancement spectra in Figure 5a,b,d are 
suggestive of the presence of multiple components in the 
intermediate-energy range. 

Discussion 

We found that both the superconducting gap anisotropy 
and the renormalized dispersion show the striking evolu- 
tion with underdoping. These behaviors are considered to 
be dependent on the extent of the screening. In associa- 
tion with the forward elastic or inelastic scatterings, the 
screening breakdown would enhance the low-energy kink. 
From the aspect of the impact of off-plane impurities, 
the inadequacy of static screening would inevitably lead 
to the nanoscale inhomogeneities, as observed by scan- 
ning tunneling microscopy experiments [34]. The forward 
scatterings by the remaining potential would generate 
additional incoherent pair excitations, as expected from 
the nodal gap suppression at low temperatures [8,25]. 



From the aspect of the electron-phonon coupling, the 
inhomogeneous depletion of the electrons for screening 
may considerably increase the coupling strength, provid- 
ing an account for the unexpectedly strong dispersion kink 
[35] and a candidate for the strong pairing interaction [8]. 
The former and latter aspects have negative and positive 
effects, respectively, on the superconductivity. Thus, we 
speculate that the doping dependence of Tq is eventually 
determined by the balance between these effects. 

Conclusions 

Summarizing, the evolution of a d-wsye high-Tc super- 
conducting state with hole concentration has been 
depicted on the basis of the high-resolution ARPES spec- 
tra of the quasiparticles and discussed in relation to 
the screening by electronic excitations. The divergence 
between the nodal and antinodal gaps can be interpreted 
as an effect of the incoherent pair excitations inherent 
in the strong coupling superconductivity. The low-energy 
kink, which rapidly increases with underdoping, should 
be caused by the forward elastic or inelastic scatterings, 
although it remains as an open question which scatter- 
ing is more dominant. The quantitative simulation of the 
doping-dependent effect will be helpful for resolving this 
problem. 
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